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ABSTRACT 



Aims. The amplitude and scaleheight of the Galactic gas disc density are not axisymmetric against expectations in a self-gravity 
axisymmetric disc. However, this lopsidedness can be explained in terms of intergalactic accretion flows, which produce non- 
axisymmetric pressure on the disc. This mechanism could be also responsible for the formation of a warp. 
Methods. We analytically derive the relationship between the disc density and the self-gravity and external pressure. 
Results. The same scenario of accretion as we proposed years ago to explain the formation of the warp explains the azimuthal depen- 
dence of the density and its scaleheight, with minimum/maximum in the positions of maximum amplitude of the warp (ip ^ 95° and 
275°), as expected from its pressure distribution. 
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c3 1. Introduction 



There are several proofs that the scaleheight of the gas disc in our 
Galaxy depends both on the galactocentric distance (Narayan & 
Jog 2002; Nakanishi & Sofue 2003) and the azimuth (Voskes 
1999; Levine et al. 2006; Kalberla et al. 2007). And there are ex- 
planations for the radial dependence of the scaleheight of galac- 
tic discs in spiral galaxies (flares) in terms of self-gravitating 
discs (e.g., Narayan & Jog 2002) or magnetic fields (Battaner 
& Florido 1995). However, it is not easy to understand the az- 
imuthal dependence in the same terms, except for an exotic pro- 
posal of the existence of a ring of dark matter embedded in the 
galactic disc with a radius that depends on the azimuth (Kalberla 
et al. 2007). We think we have a better solution that is less ex- 
otic, less ad hoc, and that agrees with other observed galactic 
features, too. 

In Ldpez-Corredoira et al. (2002, hereafter LBB), we pro- 
posed a mechanism to explain the formation of warps (S-shaped 
or U-shaped or a combination of both) in spiral galaxies: the 
accretion of the intergalactic medium (IGM) onto the disc (see 
LBB, Fig. 6). Indeed, up to now there have been no alternative 
explanations for the formation of U-warps. No massive halo is 
necessary, or high values of magnetic fields, or satellite com- 
panions are necessary, although the presence of these elements 
would not modify qualitatively the present conclusions. In a 
Milky-Way-like galaxy, the mean density of baryonic matter in 
the IGM needed to produce the observed warp is around 10"^'' 
kg/m^ when the infall velocity at a large distance is ~ 100 
km/s (LBB). These numbers were corroborated independently 
by Sanchez-Salcedo (2006). This hypothetical low-density net 
flow is a very reasonable physical assumption that would ex- 
plain why most spiral galaxies are warped. There are also dif- 
ferent types of observations successfully explained by the LBB 
hypothesis (Lopez-Corredoira et al. 2008): 1) accretion of ~ 1 



Mo/yr of low metallicity gas onto the disc as expected from the 
chemical evolution of the Milky Way (see LBB, §3.1); 2) fre- 
quency of warps and its amplitude depending on environment; 
3) lower frequency of U-warps over S-warps. 

In this paper, we claim that the same mechanism proposed 
by LBB to explain the Galactic S-warp is able to explain another 
observational fact: the azimuthal dependence of the gas density 
and scaleheight in the outer disc. 



2. Density distribution in the vertical direction of a 
self-gravitating disc 

As has been known for a long time (Spitzer 1942; Narayan & 
Jog 2002 and references therein), the distribution of matter in 
the vertical direction of the disc can be derived from applying 
the hydrostatic equilibrium. 
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and the Poisson equation. 
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where p is the mass density; (v^) the variance of the vertical ve- 
locities along the vertical z direction in an isothermal disc; 
fl0, and fl- are the three components of the acceleration per unit 
mass in cylindrical coordinates {R, 4>, z; the vertical direction 
z being defined as perpendicular to the ring); and a^^grav is the 
part of the x-component due to the self-gravity of the disc. On 
the right-hand side of Eq. (|2]), the third term is zero for an ax- 
isymmetric disc, and the second term is neglected for a thin-disc 
approximation, which are the normal assumptions in most anal- 
yses of the scaleheight. Here we do not neglect the second term 
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since the scaleheight is not small. Instead, we adopt a monopolar 
approximation for the radial acceleration. 
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gal 



«2 



(3) 



where Mgai is the mass of the Galaxy within radius R, which we 
consider constant for the outer disc. For the azimuthal acceler- 
ation a^^grav, there is some contribution from the warp, but we 
did check that it is negligible. We integrated over all the rings 
of the warped galaxy with a height ZwiR, <P), which we take from 
Levine et al. (2006), and we found that the azimuthal acceler- 
ations are of the order of 10 m/s^, and the third term of the 
right-hand side Eq. (|2]) is more than 50 times lower than the sec- 
ond term for all values of R > 15 kpc and (p. Therefore we can 
neglect this direct effect of the warp. The effect of the bar is also 
negligible: azimuthal accelerations lower than 5 x 10"'"* m/s^ 
with the potential of the long bar of Lopez-Corredoira (2007, 
Eq. 7). Therefore, if one only considers the self-gravity accel- 
eration, by doing the derivative with respect to z of Eq. ([1]) and 
assuming that the gravity is the only acceleration dependent on 
z, with Eqs. (|2| and ((S), we get 



p"-^+ Kip^ + K2P = 0, 
P 



(4) 



where Ki - K2 - The prime (') stands for the deriva- 

tive with respect to z- 

2.1. Solutions of the differential equation 

The differential equation (|4|i may be expressed as an integral. 
If we set as a boundary condition that the maximum density is 
reached at z = zo (p(zo) - A, p'(zo) - 0), we obtain 
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where x - which is a symmetrical distribution of p(z) with 
respect to zo- The half-width-half-maximum (HWHM) is 
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and the surface density is 
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For very thin discs (p very high), the second term on the 
right-hand side of Eq. (|2]i is negligible {K2 ~ 0). The solution 
to this differential equation with = is the classical squared 
hyperbolic secant solution (Spitzer 1942): piz) - A sech\a{z - 
Zo)], a = yfAKUl = KktIA. 

Since Ki and K2 do not depend on 0, it is clear that the width 
of the disc, HWHM, is independent of 0, too, and only depends 
on the radius, so the self-gravity application does not explain the 
azimuthal dependence of the scaleheight. 



3. With external pressure by a continuous accretion 
of IGIVI 

We propose that the explanation for the variation in the gas disc 
thickness depending on (p is that the external pressure due to 



the accretion of IGM onto the disc depends on 0, because the 
average flow falls down to the disc with an angle different to 71 12 
in general with respect to the plane (LBB). The pressure exerted 
over the disc would be similar to a piston mechanism, only from 
one side of the disc (Sanchez-Salcedo 2006, §4.6). The action 
is perhaps ram-pressure, due to the friction of clouds against the 
interstellar medium (Sofue & Wakamatsu 1993). 

The vertical acceleration due to this pressure, fl7,p,ei, will 
depend on z. There is a gradient of force in the vertical direc- 
tion due the higher absorption of linear momentum by the first 
layers of the disc that collide with the accreted gas. For a very 
low dense, disc, as is the case of the very outer disc, the ab- 
sorption of momentum is not total, because part of the gas can 
cross the disc completely and escape from it. We suppose that 
the external pressure is attenuated exponentially along the z-axis 

Pz,extiz) - P;,f..vf(±'^)e^'''^-'±"''*-~''*^, where C is a cross-section per 
unit mass characteristic of the interaction between the accreted 
flow and the gas disc along the vertical axis. Here we do not take 
into account the effect of the galactic rotation and the variations 
in the attenuation with different incident angles. The sign -h or - 
depends on whether the flow comes from z = c» or z = -00. The 
force per unit volume is -VP,{z). Hence, 



aW/?,0,z) = -(±)W0)C e 



±c[/;„P(zK-] 



(8) 



where F, = P, ext{±°°) is the vertical component of the external 
force per unit surface (dS - RdRdcp) due to accretion. 

We do the derivative of Eq. taking the exponential close 
to one in a low-density disc in the very outer disc [|Ccr small 
compared to one; indeed, in ^ we will see that ^Ccr <~ 
10"'^[pi,(kg/m^)]"'^^, which is smaller than one for the expected 
values of pi, {~ 10"^^ kg/rn^) although not much smaller; but we 
take this as a rough approximation]: 



-F,C^p. 



(9) 



Given a continuous inflow of particles with density at infi- 
nite distance pi,, velocity vo, and angle with respect to the plane 
60 with azimuth 0o of the direction of the inflow (so the flow 
comes from -^o, (po + n), and Galactic mass within R of Mg„i, 
the following results hold [from LBB, Eq. (45); applying Eqs. 
(28), (33), (34), (38) of the same paperfl 



F.iR, 4>)dS = 



PhvldS 



Xl 



7^2 1 sin 00 1 



(l-^-"'(^o-0)(-^-l)) 



(10) 



Rx2 



i 



R^x^ 



2 ^ RGMgaiX3 



X2 

2 



r4 gm,,,x3 

4 ^ 2vl 



R^xj RCMg„iX3 
~ + 7„ 



xi((f>) = ^J 1 - cos^ 00 sin^(0o - <P) 

X2((f>) = -\l 1 - C0S2 00 COs2(0o - 0) 
X^(<p) - 1 + COS 00 COS{<p() - <p). 

Figure [Tlplots Fj' with the parameters used by LBB: Mgai = 
2 X 10" Mq, Vo - 100 km/s, and 0o is a free parameter. Fr is 



' We have found an erratum in LBB: Eq. (28) should have opposite 
sign, i.e. eoQ = - cos(vo, rg) = - cos(0o) cos(0o - <(>)■ However, numbers 
and orientation of the warp are correct as stated in LBB. 
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Fig. 1. ' according to eq. (fTOl i with the parameters given by 
LBB andVo = 275°. 

~ 10""^ kg/m/s^. We take 0o = 275°, Oq < 0, the corresponding 
direction of the accreted flow, in order to produce the Galactic 
S-warp with U-warp southwards (LBB) and maximum of the 
S-warp at = 95°, 275°, as observed (Voskes 1999). The de- 
pendence is plotted for R - 25 kpc and variable (p, while other 
values of R give different amplitudes in the variation, but sim- 
ilar azimuthal dependence. The amplitude of the variation with 
the azimuth is strongly dependent on 9^). Clearly, for any value 
of 9q, the minimum pressure is for (p - (pQ + n - 95°, and the 
maximum is for (p - (f>o - 275°. Nonetheless, if we took into 
account the rotation of the galaxy and an external pressure atten- 
uation dependent on this and the incident angle, or the degree of 
dumpiness of IGM and disc, the diff'erence in the position of the 
maximum and minimum would not be strictly tt, and the shape 
of Fig. [T]would vary. 

3.1. Differential equation 

If we now consider that the total acceleration is a- — a^^grav + 
<^z,pres in Eq. ([TJ, assuming an immediate response of the disc to 
the pressure exerted by the accreted material (Ldpez-Corredoira 
et al. 2008), together with Eqs. ^ and (|9]l, we again get the same 

differential equation dUi but with constants K\ - ^^^^t^7^4^. 



K2 = 



CM,,, 



The total acceleration is null at zo, where 



the pressure acceleration is compensated by other gravitational 
forces to keep the warp static. With C ~ 10^ - 10^ m^/kg, we 
get values of F.{R, <p)C^ of the order of A-ttG, that is, an external 
pressure as significant as the self-gravity. And, most important, 
HWHM will depend on because depends on (p and also on 
the maximum amplitude of the density. A, and the dispersion of 
velocities depends on (f> through its dependence on the pressure 
at z = zo P(<py- 



A((P) = MM 



{vim = <y?)(0o) 

P(0) oc (l + 



f(0) 



(11) 
(12) 
(13) 



This last proportionality of the pressure stems from the exter- 
nal pressure acting like an extra self-gravity, and the pressure 
is proportional to the total acceleration. For a monoatomic gas. 



Fig. 2. Dependence of HWHM on cr, from Eqs. (|6|l and (|2l), for 
different values of Ki with the minimum value K\x) - ^ given 

when only gravity is present and other values when an external 
pressure is added, and fixed K2 and (yj). 



7 = 5/3 for adiabatic compression, and y - I for isothermal 
one. This applies in the distribution of pressure as a function of 
0; for the vertical dependence, we have already assumed that it 
is isothermal. 

In Fig. |2] we see how the HWHM is reduced for a given cr 
when K\ is increased, i.e. when F, is increased. The parameters 
are = 25 kpc, M^„/ = 2x10" Mq, and cr,,^ = 10 km/s. 



4. Comparison with the observations 

Voskes (1999, Fig. 15) and Levine et al. (2006, Fig. 5) have 
shown that the scaleheight of the outer disc {R > 20 kpc) is 
2-3 times higher on average for < <p < 180° than for 
180° < (p < 360°. A more accurate estimation of the maximum 
of the scaleheight is derived by Kalberla et al. (2007, Figs. 18, 
19), who place it for 90° < (p < 110°, while the minimum for 
250° < 4> < 270°. Within the LBB scenario, two possible di- 
rections are possible to produce the observed S-warp, and this is 
one of them: the wind coming from the direction of the northern 
warp, although this solution could not explain the asymmetry 
of the southern/northern warp as a sum of S-hU warp. A higher 
pressure is expected for the region around the southern warp and 
consequently a lower thickness therein. Not only is the pressure 
lower at ^ 90°, but the surface density cr is also lower than for 
the average value of criR) (Voskes 1999, Fig. 13; Levine et al. 
2006, Fig. 1; Kalberla & Dedes 2008, Fig. 9) and the amplitude 
of the density A (Kalberla & Dedes 2008, Fig. 8). This is another 
fact explained by our model by means of Eqs. Q and ( fTTT i: the 
lower the pressure the lower the density. 

If we take the values observed by Kalberla et al. (2007, 2008) 
at = 25 Kpc, excluding the bin of 90° < < 110°, which 
is possibly an outlier in a special region: Aijpo) = 1.5 x lO"^'* 



kg/m 



3. A(0o+a-) 



7, 



a"(0o) 



HWHM{<l>a) 



2.5. This 



A((Ao) ' o-(0o) ' HWHM{<Po+n) 

last number is not independent of the other three. The param- 
eters that better fit these numbers are: y - 1.22, 9(, - -31°, 
PhC^ - 1.8 X 10 m/kg, giving an azimuthal dependence of 
HWHM, A and cr as plotted in Fig. |3] The average accretion at 
R^25 kpc is 5 X lO"'' [w^H^) kg/m^/Gyi". which is around 2 
times the average surface density per Gyr; note, however, that 
only a small ratio of its linear momentum (~ Ccr) is transmitted 
to the Galactic disc. 

For the shape of the scaleheight as a function of for 
R - 25 kpc (or any other value of the galactocentric distance 
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Fig. 3. Azimuthal dependence of the variables HWHM, A and cr 
for R - 25 kpc. Parameters: vo - 100 km/s, Mgai = 2 x 10" 
Mo, cr,^{<po) = 10 km/s, A(0o) = 1.5 x lO^^^ kg/m\ 0,) = 275°; 
r = 1.22, 00 = -31°, piC^ = 1.8 X lO-'"* m/kg. 
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20 < < 40 kpc), we observe in Fig. 19 of Kalberla et al. (2007) 
a strongest dependence of the scaleheight with (p near the max- 
imum; however, the variation in amplitude is very slight around 
the minimum of flaring. This is not observed in the shape of Fig. 

El 

Therefore, at present with our simple model, we cannot ex- 
plain the exact azimuthal dependence. Possibly, our rough as- 
sumptions in the present and LBB calculations need to be im- 
proved to get a more accurate result, assuming that Kalberla 
et al. (2007, 2008) are also correct with their assumptions on 
the kinematics of the Galaxy. Indeed, the Kalberla et al. (2007, 
2008) density distributions are derived by assuming a non- 
axisymmetric disc using epicyclic corrections in the rotation 
curve, consistent with their hypothesis to explain the asymme- 
tries in the outer disc. Since we can explain the asymmetries of 
the outer disc with the accretion of IGM flows, it is possible that 
those epicyclic corrections are not necessary, although they are 
not necessarily inconsistent with our scenario. In any case, our 
model describes the asymmetries roughly. The value of the pa- 
rameters for the infall material also more or less agree the values 
necessary for the warp production in LBB. If Kalberla et al. as- 
sumed axisymmetric rotation, they would get a gas distribution 
that is several kpc more extended to the south (180° < / < 360°) 
than to the north (Levine et al. 2006, Fig. A18). This would also 
agree with our scenario because we predict a higher density in 
the south in general. 



5. Conclusions 

The same LBB scenario that explains the formation of warps is 
able to directly explain the azimuthal dependence of the HI den- 
sity distribution in the Milky Way, with maximum and minimum 
in the positions expected, without introducing ad hoc elements. 
Sanchez-Salcedo (2006) raises the criticism that LBB mecha- 
nism is not plausible because it would produce a dependence on 
the scaleheight of the disc with the Galactocentric azimuth, 0, in 
the outer disc. Rather than being an objection, it is another argu- 
ment in favour of our model because this dependence is actually 
observed in our Galaxy. 
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